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ABSTRACT 



Context. Long-duration gamma-ray bursts (LGRBs), which release enormous amounts of energy into the interstellar medium, occur 
in galaxies of generally low metallicity. For a better understanding of this phenomenon, detailed observations of the specific properties 
of the host galaxies (HG) and the environment near the LGRBs are mandatory. 

Aims. We aim at a spectroscopic analysis of HG 031203, the host galaxy of a LRGB burst, to obtain its properties. Our results will be 
compared with those of previous studies and the properties of a sample of luminous compact emission-line galaxies (LCGs) selected 
from SDSS DR7. 

Methods. Based on VLT/X-shooter spectroscopic observations taken from commissioning mode in the wavelength range ~AA3200 
- 24000A, we use standard direct methods to evaluate physical conditions and element abundances. The resolving power of the 
instrument also allowed us to trace the kinematics of the ionised gas. Furthermore, we use X-shooter data together with Spkzer 
observations in the mid-infrared range for testing hidden star formation. 

Results. We derive an interstellar oxygen abundance of 12 -I- log O/H = 8.20 ± 0.03 for HG 031203. The observed fluxes of hydrogen 
lines correspond to the theoretical recombination values after correction for extinction with a single value C(Hj8) = 1.67. We produce 
the CLOUDY photoionisation H n region model that reproduces observed emission-line fluxes of different ions in the optical range. 
This model also predicts emission-line fluxes in the near-infrared (NIR) and mid-infrared (MIR) ranges that agree well with the 
observed ones. This implies that the star-forming region observed in the optical range is the only source of ionisation and there is no 
additional source of ionisation seen in the NIR and MIR ranges that is hidden in the optical range. We find the composite kinematic 
structure from profiles of the strong emission lines by decomposing them into two Gaussian narrow and broad components. These 
components conespond to two H ii regions, separated by ~34 km s"', and have full widths at half maximum (FWHM) ~1 15 and ~270 
km s"', respectively. We find that the heavy element abundances, extinction-corrected Ha luminosity L(Hq')=7.27 x lO'" erg s"', 
stellar mass M, = 2.5xlO*Mo, star-formation rate SFR(Hq') = 5.74 Mq yr"' and specific star-formation rate SSFR(Hq') = 2.3x10"^ 
yr"' of HG 031203 are in the range that is covered by the LCGs. This implies that the LCGs with extreme star-formation that also 
comprise green pea galaxies as a subclass may harbour GRBs. 
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^ . 1. Introduction 



GRB 031203 is one of the two closest (z < 0.2) long-duration 
gamma-ray bursts (LGRB) known apart from the exceptional 
GRB 980425. 

Most soft-spectrum LGRBs are accompanied by massive 
stellar explosions (the GRB-supernovae (SNe) connection, see 
van Paradiis', 799?: 'W ooslev & Blooml l2006t ISoderbergl l2006t 
Hiorth & Bloom, 2011). This and other evidence links LGRBs 
to ongoing star formation and consequently to the most massive 
stars as possible progenitors of LGRB s and X-ray flares (XRFs) 



conclude that there is no compelling evidence that HGs are 
peculiar galaxies and that they are instead similar to nor- 
mal star-forming galaxies in the local and distant universe. 
Accordingly, HGs are proposed to be used as tracers of star 
formation. Spectroscopic studies of HGs give information on 
LGRB progenitors and on the physical properties of regions 
hosting LGRBs . Investigations of HGs have been performed for 
many galaxi es (Christensen et al.L 20041 : Prochaska et al., 2004 
Stanek et"al1 , i2006.; .Sollerman et al.L 2005c iWiersema et al 



20071: iMargutti et aU 120071: iThone et all llOm iHammer et al' 



(iBlooni et al.L l2002t iLe Floc'h et al.L I2003L IChristensen etal. 

I2004L iFruchter et al.L l2006h . Wolf-Ravet (WR) spectral fea 

tures w ere found in th e spectra of sever a l LGRB host galaxies iSchadv et al.L 201 
(HGs) (iHammer et al 1 1200(a iHan et al.L l2010l) . LSavagUo et all 
(l2009l) present an extensive study of the most extensive sam- 
ple of HGs to date, encompassing 46 targets. The authors 



2006; Kewlev et al.','2007^,'Th6ne et al.V2008'; 'Christensen et al 
2008; Savaglio et al., 2009; Han et al. ,2010; Syen sso n et al 
20101: iWatsonetaU |20 1 OL iLevesgue et al.L 120101 



2011 
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* Based on observations collected at the European Southern 
Observatory, Chile, ESO programme 60.A-9024(A). 



IVerganietaLL 1201 iL IWiersemaL l201ll) 
However, to date spatially resolved stu dies have been per- 
formed only for severa l low-redshift H Gs (Hammer et al.', "20061 
Thon e etal.. 2008; Christensen et all 2008; Sollerman et al], 
2005L ILevesgue et al.L |20T l]). For more distant emission-line 
galaxies at redshifts z ^ 0.1, the properties of the LGRB environ- 
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Fig. 1. X-shooter acquisition image of HG 031203 [ESO pro- 
gram 60.A-9442(A)] showing the slit location. The galaxy has a 
compact, slightly elongated shape and likely consists of two H n 
regions. North is up and East is to the left. 



ment can be retrieved only from the interstellar medium (ISM) 
of the entire galaxy, as is the case for HG 031203. 
(1201 Oh. 



Svensson et al 



evesQue et al 



(I2010ah 



ISavadio et al.| (l2009h . fKewlev et al.l (120071) and I Wiersema et al.l 
( l2007h found that HGs are relatively low-mass galaxies 
with a low-metallicity ISM and intense star formation. For 
a given mass (or luminosity) HGs are systematically offset 
towards lower m etallic iti es in lumino si ty-met al Ucity (L - Z) 



(iMargutti et al.L 120071: 



(IMargu 
I2007t I: 



iLevesQue et al 



('Savagli o et all 12009 : 



r2010bt IVergani et al 



Stanek et all l2006t iKewlev et al 
2010ah orlna s s-meta l licity (M - Z) 



■Svensson et al.L 120101: iLevesque et al . 
I20T ij) relations as compared to dwarf 



1995 



2002; 



irregulars and normal star-forming emission-line galaxies 
dLequeux et al.L [19791 I Skillman et al.', 198?; 'Richer & McCalll 

Kobulnicky & Z aritskv. 1999; Melbourne & Salzen, 

iLeeetaU 12001 IPi vugin et all 120041: iKongL l200l 



Shi et all 12001 iLee et al.L l2006l) . Moreover, ILevesque et al.l 
(I2010bh obtained the offset between nearby HGs with z < 0.3 
and nearby SDSS star-forming galaxies as well as between 
HGs with intermediate-redshift (0.3 < z < 1) and emission- 
line galaxies from the DEEP2 survey (<z>=0.8). It was also 
established that LGRB HGs have lower metallicities t h an SN 
HGs withou t accompanying LGRBs ( Fruchter et all 20061: 
[Modiaz et al.', 2008; Levesaue et al.L l2010al) . iKobulnic kv et alJ 
42003) and Savaglio et al. (2005i) have found a temporal evolu- 
tion of emission-line galaxies in L - Z and M - Z diagrams, 
respectively. Higher redshift galaxies are more metal-poor at 
a fixed mass or luminosity than local emission-line galaxies. 
Host galaxies (HGs) are characterised by high star-formation 
rates (SFRs) a n d hig h specific star-formatio n rates (SSFRs) 
( Savaglio et al.L 120091: ISvensson et al.L 120101; ILevesque et al 



I2010aibh . 

The metallicities of most HGs are obtained using empiri- 
cal strong-line methods, such as R23, 03N2, rn a ss-me t allicity 
relation and o t her ca lib rations ([Savaglio etal 1_|2005", 12001 
Levesaue et al.L l2010allH : ISvensson et al.L |2010L F^rgani et al^, 
201 ll) and are somewhat uncertain. This is because of a well- 



known offset between the metallicity obtained by the direct Tg- 
method and empirical strong-line methods. Oxygen abundances 
obtained by empirical methods are by ~ 0.2-0.6 dex higher 
than those obtained with th e -m ethod dGuseva et al.L I2009L 
iHovos et all I2005L IShi et all l2005l) . Moreover, the luminosity- 
metallicity relation is also blurred because the metallicity de- 
terminations of various galaxy samples do not employ a unique 
technique. Therefore, we use only a comparison sample of HGs 
for which the metallicity is derived using the r4.-method. 

The metallicity of the ISM in the region around the LGRB 
site may be different from that averaged over the whole galaxy. 
However, oxygen abundance variations over dwarf compact 
galaxies are small (~0.1 dex), that is, in the range of errors 
of abundance determ inations (see, e.g. iPapaderos et al.L I2006L 
llzotov et"an. l2006bh . Therefore, the average metallicity of the 
dwarf HG is comparable to that of the LGRB environment. Thus, 
in low-mass galaxies the metallicity of an entire galaxy can be a 
proxy of an LGRB environment. 

Therefore, a comprehensive study of HGs that exhibit bright 
H II region features is essential for a better understanding of 
LGRBs and needs to investigate the physical conditions, chem- 
ical abundances, extinction, and kinematic structure of their en- 
vironment. 

The long-duration (~ 30 s) GRB 03120 3 was discovered in 
2003 with INTEGRAL (iGotz et all l2003l) . A compac t dwarf 
galaxy coinciding with the X-ray source (iHsia et al 1 I2OOI 
was later id entified as the GRB h ost HG 031203 at redshift 
z - 0.1055 dProchaska et al.L |2004 . Monitoring of the galaxy 
to search for SNe eve nts led to the discovery of SN20031w 
(iTagliaferri et al.L [20041) . The proximity of the galaxy (it is one 
of the closest known long-duration GRB host galaxies) allows 
us to perform a detailed analysis of its properties. 

An oxygen abundance 12 + log O/H = 7.96 (r^-method) 
and high reddening E{B - V) - 1.17 co r respond ing to C(Hj6) 
- 1 .72 were obtained by ILevesque et al.l ( l2010a ) from Keck I 
LRIS observations. Based on optical VLT spectra Margutti et al.l 



([2OOI derived an oxygen abu ndance 12 + log O/H = 8.12 
+ 0.04. iProchaska et"alT(l2004l) from Magellan/IMACS optical 
spectroscopy obtained an oxygen abundance 12 + log O/H = 
8.02 + 0.15. This galaxy is located at low galactic latitude, 
therefore the extinction by the Milky Way is high. Special at- 
tention is needed to derive background and intrinsic extinc- 
tion in HG 031203, which c an a ffect the abundance de termi- 
nation. iMargutti et al.l (l2007h and IProchaska et all (|2004 anal- 
ysed both the Milky Way and the intrinsic extinction in HG 
031203 and derived intrinsic extinction coefficients C(Hy6) = 
0.59 and C(H/?) = 0.46, respectively. Their intrinsic host-galaxy 
extinctions seem higher that usually obtained in star-forming 
compact dwarf galaxies. This is mainly because the authors as- 
sumed a lower value of the Galactic reddening than that given by 
ISchlegel et al.l (Il998h . On the other hand, the available literature 
data show no evidence f or substantial ext i nction in any compre- 
hensi vely studied HGs dSoUerman et al.L I2005L ISavagfio et all 
2009). Nine o ut of e leven well-studied GRB HGs collected by 
Savaglio et al., d2009l) have A{V) in the range -0-0.6 [or C(Hy8) 
~0-0.3] as deri ved from the B aimer de crement. For HG 030329 
and HG 980425 lSollerman et al.l d2005l) derived C (H6)Hr, of 0.06 
and 0. 10, re spectively. Therefore, i t is likely that IMargutti et al.] 
(l2007h and IProchaska et al.l d2004l) overestimated the intrinsic 
reddening in HG 031203. 

We here present high-quality archival VLT/X-shooter spec- 
troscopic observations of HG 031203 over a wide wavelength 
range, which allows us to derive physical conditions, intrinsic 
reddening and element abundances in the ionised interstellar 
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Fig. 2. Flux-calibrated VLT/X-shooter near-UV and optical range (UVB + VIS arms) spectrum of the GRB 031203 host galaxy 
corrected for a redshift of z = 0.1055 (upper spectrum in each panel). The lower spectrum is the upper spectrum downscaled by a 
factor of 100. The scale of the ordinate is that of the upper spectrum. 



m edium of the galaxy . The same observations were discussed 
bv lWatson et al.lOOlOh . However, they did not draw any conclu- 
sions on extinction and heavy element abundances. The observa- 
tions encompass the near-infrared (NIR) range that allows us to 
find out whether hidden star formation is present in this galaxy. 
Such a broad wavelength range extending to the NIR range per- 
mits us to derive a more reliable stellar mass (and consequently 
SSFR) of the galaxy. This is because cool low-mass stars mainly 
contribute to the stellar mass and they also emit mainly in the 
NIR range. We note, however, that in low-metallicity compact 
dwarf galaxies with high SFR the presence of strong ionised 
gaseous emission complicates the stel lar mass determinat ion. 
We also use the Spitzer observations of IWatson et al.l (l20I0h for 
testing more hidden star formation in the mid-infrared (MIR) 
range. 

The observations are described in Sect. |2] In Sect. [3] we 
present the main properties of the host galaxy. In particular, the 
diagnostic diagram is given in subsect. 13.11 The element abun- 
dances in HG 031203 are presented in subsect. [32] The extinc- 
tion and hidden star formation are discussed in subsect. 13.31 
the CLOUDY modelling and the comparison of the observed 
and predicted emission-line fluxes in the MIR range are dis- 
cussed in subsect. 13.41 the kinematic structure in subsect. 13.51 
the luminosity-metallicity relation in subsect. 13.61 and the star- 
formation rate in subsect. 13.71 Our conclusions are summarised 
in Sect.S] 



2. Observations 

A new spectrum of the GRB 031203 host galaxy was obtained 
during the commissioning of the VLT X-shooter on 2009 March 
17 [ESO program 60.A-9024(A)]. In Fig. [T]the acquisition im- 
age of HG 031203 with the slit location is shown. The galaxy 
presents itself as a slightly elongated source with major axis 
~ 1", which is comparable to the slit width of 0'.'9 - 1". Two 
H II regions are present in the galaxy, one brighter than the 
other. The observations were performed in the wavelength range 
~/l3200 - 24000A at the airmass 1.037 using the three-arm 
echelle X-shooter spectrograph mounted at the UT2 Cassegrain 
focus. In the UVB, VIS and NIR arms the total exposure time of 
4800s were broken into four equal subexposures of 1200s each. 
Nodding along the slit was performed according to the scheme 
ABBA with the object positions A or B differing by 5" along the 
slit. In the UVB arm with wavelength range /13233 - 5600A a slit 
of l"x 11" was used. In the VIS and NIR arms with wavelength 
ranges ^5475 - 10206A and ill 100 - 23985A, respectively, a 
slit of 0'.'9 X 1 1 " were used. The binning factor along the spa- 
tial and dispersion axes in the UVB and VIS arms was 1, except 
for the UVB arm where the binning factor along the dispersion 
axis was 2. The above instrumental set-up resulted in resolving 
powers A/AA of 10200, 8800, and 5100 for the UVB, VIS, and 
NIR arms, respectively. The seeing was ~ 1". All observations 
were obtained at low airmass of ~ 1 .03 - 1 .06, therefore the ef- 
fect of the atmospheric dispersion is low. The correction for this 
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Fig. 3. Same as Fig.|2]but for the near-infrared (NIR) range. 




2.05 



effect was applied to the UBV and VIS arm spectra during ob- 
servations. We used the spectrum of the spectrophotometric stan- 
dard star GD 153 for the flux calibration. The absolute synthetic 
fluxes for this star were taken from the Space Telescope Science 
Institute web site Q. Spectra of thorium-argon (Th-Ar) compar- 
ison arcs were used for wavelength calibration of the UVB and 
VIS arm observations. For the wavelength calibration of the NIR 
spectrum we used night-sky emission lines. 

In a first step cosmic ray hits of all UVB, VIS, and NIR 
spectra were removed using the routine CRMEDIAN. The re- 
maining hits were subsequently removed manually after back- 
ground subtraction. The two-dimensional UVB and VIS spec- 
tra were bias-subtracted and flat-field corrected using IRAlQ. 
The two-dimensional NIR spectra were corrected for dark cur- 
rent and divided by the flat frames to correct for the pixel sen- 
sitivity variations. For each of the UVB, VIS, and NIR arms 
we separately coadded spectra with the object at the position 
A and spectra with the object at the position B. Then, the coad- 
ded spectrum at the position B was subtracted from the coad- 
ded spectrum at the position A. This resulted in a frame with 
subtracted background. We used the IRAF software routines 
IDENTIFY, REIDENTIFY, FITCOORD, and TRANSFORM to 



' ftp://ftp.stsci.edu/cdbs/current_ calspec/ 

- IRAF is the Image Reduction and Analysis Facility distributed by 
the National Optical Astronomy Observatory, which is operated by the 
Association of Universities for Research in Astronomy (AURA) under 
cooperative agreement with the National Science Foundation (NSF). 



perform wavelength calibration and correct for distortion and 
tilt for each frame. The one-dimensional wavelength-calibrated 
spectra were then extracted from the two-dimensional frames us- 
ing the APALL routine. We adopted extraction apertures of l"x 
2'.'5, a:9 X 2/5 and a/9 x 2/5 for the UVB, VIS and NIR spec- 
tra, respectively. Before extraction, the spectra at the positions 
A and B in the two-dimensional background-subtracted frames 
were carefully aligned with the routine ROTATE and co-added. 
For the flux calibration we adopted an atmospheric extinction 
curve for the Cerro Tolol o Observatory, wh ich is very similar 
to that for Cerro Paranal (iPatat et al.L l20Ilh . Although telluric 
standards were also observed, no attempts were made to correct 
for telluric absorption. This is because this correction resulted in 
very noisy spectra in the regions of strong absorption and does 
not help much in the determination of the reliable fluxes of the 
lines that fall in those wavelength ranges. Fortunately, there are 
many strong lines in regions of low or no telluric absorption. 
The number of these lines is sufficient for our analysis. In par- 
ticular, the redshifted strongest NIR emission line Pao- is in the 
region without telluric absorption. The resulting flux-calibrated 
continuum in all orders is monotonic over the entire range of the 
spectrum. Therefore, no adjustment of the UVB, VIS, and NIR 
spectra were needed. 

The resulting flux-calibrated and redshift-corrected UVB 
and VIS spectra of HG 031203 are shown in Fig. |2] and the 
resulting flux-calibrated and redshift-corrected NIR spectrum is 
shown in Fig. [3] 
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Fig. 4. Baldwin-P hillips-Terlevich (BPT) diagram 
([Baldwin et al.L |1981[) for luminous com pact galaxies (LCGs) 
from SDSS DR7 (llzotov et al.L 1201 lab (small blue circles). 
HG 031203 (this paper) is shown by a large filled r e d star . 
The five low-me t alficity AGNs from llzotov & Thuaril (l2008h 
and llzotov et"al] (1201 Oh are indicated by large filled black 
circles. Host galaxy (HG) 031203 for any data is shown by a 
star Other HG s are d enoted as follows (see t ext for details): 
ISavagUo et af] ( l2009l) : filled green squares. iLevesaue et al.l 
(I2010ah : filled light blue squares, (IChristensen et all l2008l) : 
violet asterisks. Hammer et al. (2006J: large green open circles 



dWiersema et al Ll2(307h: fifled purp l e trian gle! iHan et al.l (1201 Oh : 
purple diamond s, Margutti et al.l (l2007l): filled purp le star. 



iProchaska et al.l ( |2004|) : open bltie star, IWatson et al. ( 1201 Ol) : 
open yellow star. All HGs with emission line fluxes lower than 
2 X 10"'^ erg s"' cm"- are shown by small symbols. Also, 
the 100,000 emission -line galaxies from SDSS DR7 selected by 
llzotov et"al] (1201 lab are shown b y grey dots. The dashed line 
represents the empirical line of iKauff'mann et al.l (l2003l) that 
separates galaxies dominated either by star formation or by an 
AGN. The continuous line shows t he upper limit f o r pure ly 
star-forming galaxies (SFGs) from IStasiiiska et alJ (l2006l) . 

(A colour version of this figure is available in the online 
journal.) 

3. Results 

3.1. Diagnostic diagram 

The position of HG 031203 in the [O iii] ^5007/^8 vs. 
FN n1 /16583/H g Baldwin-Phillips-Terlevich (BPT) diagram 
tealdwinetaU 11981!) is shown in Fig. |4]by a large filled red 
star. Lumin ous compact emission-line galaxies (LCGs) from 
SDSS DR7 (llzotov et al.L 1201 lah are shown by small blue cir- 
cles. The global properties of LCGs are similar to those of 
the green pea galaxies, which were recently discovered by 
[Cardamone et al. (2009) as a new class of luminous compact 
galaxies at redshift z ~ 0.1 - 0.3 with a peculiar bright green 
colour. Based on the photometri cally selected sample of 251 
galaxies, ICardamone et al.l (l2009l) found that green pea galax- 
ies have some of the highest specific star-forming rates seen 
in the local Universe and are similar in size, mass, luminos- 
ity, and metallicity to luminous blue compact galaxies and to 
high-redshift ultraviolet l uminous galaxies (Ly man-break galax- 
ies and Ly-ff emitters), llzotov et all ( 1201 la!) constructed and 
studied an extensive sample of 803 star-forming LCGs in a 
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Fig. 5. Abundance ratios log N/O (a), log Ne/O (b), log 
S/O (c) and log Ar/O (d) vs. oxygen abundance 12 -i- log 
O/H f or HG 0312 3 (this paper) and for LCGs from SDSS 
DR7 (llzotov et all 1201 lah are shown by large red filled 
stars and small blue filled circles, respectively. Data from 
the literature and recalculated from published emission line 
fluxes are shown by filled and open symbols, respectively. 
All data for HG 031203 are shown b y stars of different 
colours. S pecifically, Margutti et al. I dlOOTh : filled and open pur- 
ple stars. Hammer et al.l (l2006l) : filled and ope n black circles, 
Wiersema et al. (2007): filled yel low triangles, Levesaue et al.l 
(I2010ab : open light blue squar es, iHan et al" I (I2010h : open blue 
rhombs. IProchaska et al.l(l2004|): filled green stars (see text). The 
solar abundance ratios bv lAsplund et al.l (l2009l) are indicated 
by the large purple open circles and the associated error bars. 

(A colour version of this figure is available in the online jour- 
nal.) 



wider redshift range z - 0.02 - 0.63, selected from SDSS 
DR7 with the use of spectroscopic and photometric criteria. 
The fi ve lo w-metal licity AGN candidates from Izotov & Thuaiil 
(l2008l) and llzotov et all (l2010l) are shown by large filled black 
circles. Data of other authors for HG 031203 are shown by 
stars in different colo urs. Other HGs are denoted as follows: 
ISavagUo et all (l2009h : filled green squares (031203, 020903 
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and 060505). lEevesque et all (l2010ah: filled hght blue squares 
(030329, 020903 and 060218'). lchristensen et a"ill2008h : violet 
asterisks (980425 - W R, SN regions and entire host galaxy), 
iHammer et al ] (I2006h : large gre en open circles (020903 and 
980425 (WR and SN regi ons')! dWie rsema et all l2007h : filled 
purple triangle (060218), jHaneOD" (2010): purple diamonds 
(020903, 03 0329, 031203, 06021 8, 060505 GRB site and en- 
tire galaxy). iMargutti et al.1 (l2007h : filled purple s tar (031203) 
iProchaska et al.l ( l2004t) : open blue star (031203). IWatson et alj 
(120101) : open yellow star (031203). All HGs with emission line 
fluxes lower than 2 x 10"'^ erg s ' cm"^ are shown by small 
symbols. Also, the 100,000 emission-line galaxies from SDSS 
DR7 are seen as a cloud of grey dots. The dashed line repre- 
sents the empirical line of iKauffmann et ahl ( l2003h . which sep- 
arates star-forming galaxies and AGNs. The continuous line 
delineates the upper limit for pure star-forming galaxies from 
IStasihska et"al](l2006l) . 

In Fig. |4] almost all LCGs lie in the low-metallicity part of 
what is usually considered as the region of star-forming galax- 
ies. HG 031203 is situated on the border between SFG and 
AGN regions and is close to the region of very low-metallicity 
AGNs. Photoionisation models of AGN show that lowering 
their metallicity moves them to the left of BPT diagram, so 
that they end up in t he SFG region (Stasihska et al., 2006). 
lIzotov&ThuanI (12008) have concluded that the BPT diagram is 
unable to distinguish between SFGs and low- metallicity AGNs . 
The five low-me t allicity AGN candidates (llzotov & ThuanL 
l2008t llzotov et"an , 1201 Ol) also contain the broad components of 
strong emission lines Ha and HyS with high luminosities and very 
steep Balmer decrement. The position of AGN candidates in 
BPT diagram can be accounted for by the models with the non- 
ther mal radiation cont r ibutin g less than ~ 10% to the total radia- 
tion (llzotov & ThuanL l2008l) . In particular, the pure photoion- 
isation model without nonthermal radiation agre es quite well 
with o bservations (see subsect. 3.4 in this paper). IWatson et aP 
(1201 Ol) rule out the presence of an AGN in HG 031203 based on 
flux ratios of optical emission lines (X-shooter data) and mid- 
IR spectrum. The auth ors confirm the previous conclusions of 
IProchaska et al.l (l2004l) andlMargutti et al.1 (l2007h . On the other 
hand. iLevesque et al.l ( 1201 Oal) . using hne-flux ratios from the 
Keck I LRIS spectrum, concluded that HG 031203 shows evi- 
dence of AGN activity. 

Overall, HGs including HG 03 1203 as well as the five low- 
metallicity candida te AGNs from llzotov & ThuanI (l2008h and 
llzotov et al.l (1201 Oh lie in the region occupied by LCGs, imply- 
ing the similarity of their properties. The location of HG 031203 
in the upper left region of BPT diagram can be accounted for 
by the young age of the starburst. Starbursts with older age are 
instead located in the right lower part of the SFG region. 

3.2. Element abundances 

We derived element abundances from emission-Hne fluxes us- 
ing a classical semi-empirical method. These lines trace the ISM 
of HG 031203. The fluxes in all spectra were measured using 
Gaussian fitting with the IRAF SPLOT routine. The line flux er- 
rors include statistical errors derived with SPLOT from non-flux- 
cahbrated spectra, in addition to errors introduced by the abso- 
lute flux calibration, which we set to 1 % of the line fluxes, ac- 
cording to the uncertainties of absolute fluxes of relatively bright 
standard stars (IOkail990l:IColina & BohlinLll994l:lBohhnLll996t 



Given a function f{x,y, ...,z), the uncertainty cr(/) is calculated 

as 



llzotov & Thuani l2004 . These errors will be propagated below 
into the calculation of the electron temperatures, the electron 
number densities, and the ionic and total element abundances. 



cr(f) 



(g.(x))\(|.(,))\....(|.(z)^' 



(1) 



The fluxes were corrected for e xtinction, using the red- 
dening curve of ICardelli eta l. (198^ and fo r underlying hy- 
drogen stellar absorption ( Izotov et all Il994l) . The equivalent 
width of the Hy6 emission line is EW(Hj0) = 134 A correspond- 
ing to an age of 3-4 Myr At this age the predicted equiv- 
alent width of the Balmer absorption lines is EW(abs) ~ 2- 
3A (IGonzalez Delgado et all 119991: iGuseva et all l2003l) wifli 
relatively small variations between diff'erent hydrogen lines. 
Therefore we assume that EW(abs) is the same for all hydro- 
gen lines in the same galaxy but varies from galaxy to galaxy. 
This assumption is justi fied by the evolutionary stellar popula - 
tion synthesis models of IGonzalez Delgado et al.l (Il999l I2d05h . 
We also note that 20 - 30% uncertainties in the EW(abs) for the 
3-4 Myr starburst result in < 0.2% uncertainties in the fluxes 
of the strongest emission lines HyS and Ha, which are most im- 
portant for the determination of the extinction coefficient. The 
extinction coefficient C(HyS) and equivalent widths of the hy- 
drogen absorption lines EW(abs) are calculated simultaneously, 
which minimises the deviations of corrected ffuxes I{A)/I(H/3) of 
all hydrogen Balmer lines from their theoretical recombination 
values as 



F(A) EW(A) + EW(abs) 



EW(B/3) 



EW(A) 



£W(H/3) + EWiabs) 



Here f(A) is the reddening function normalised to the value at 
the wavelength of the H/3 line, F{A)/F(H/3) are the observed 
hydrogen Balmer emission line ffuxes relative to that of H/?, 
EW(/l) and EW(H/3) the equivalent widths of emission lines, 
and EW(abs) the equivalent widths of hydrogen absorption lines. 
The extinction-corrected ffuxes of emission lines other than 
those of hydrogen are derived from the equation 



1(A) 



F(A) 



lim Fim 

The derived C(Hy6) is applied for correction of all emission- 
line ffuxes in the entire wavelength range AA32QQ - 24000A. 

We varied Ry, the ratio of total to selective extinction, 
in the range from 2.5 to 4.0 and found that these variations 
change the relative line ffuxes by less than ~2%. In the follow- 
ing we adopt Ry - 3.2. The extinction-corrected relative ffuxes 
100x/(/l)//(HyS) of the lines, the extinction coefficient C(Hy6), the 
equivalent width of the H/? emission line EW(Hy6), the Hy6 ob- 
served ffux F(Hy6), and the equivalent width of the underlying 
hydrogen absorption lines EW(abs) are given in Table [1] We 
note that ffuxes of Balmer hydrogen emission lines corrected 
for extinction and underlying hydrogen absorption (column 3 
in Table [Hi are clo s e to th e theoretical recombination values of 
iHummer & StorevI (Il987l) (column 4 of the Table), suggesting 
that the extinction coefficient C(HyS) is derived reliably. 

The physical conditions and the ionic and total heavy ele- 
ment abundances in the ISM of the GRB 31203 host galaxy 
were derived following llzotov et al.1 (l2006ah (Table |2]i. In par- 
ticular, we adopt the temperature Tf''{0 iii) for 0~^, Ne^^, and 
Ar^+, which is directly derived from the [O iii] /14363/(/l4959 + 
/15007) emission-line ratio. The electron temperatures T^iO ii) 
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Fig. 6. Left panel. Decomposition of strong emission-line profiles into two Gaussian components in the spectrum of HG 03 1203 for: 
a) HyS /14861; b) [O iii] ^4959; c) [O m] ^5007; d) Yia X 6563; e) [S iii] i9532 and f) Pao- ^18756. The observed spectrum and the 
fit are shown by black solid and red dashed lines. The two Gaussian components and residual spectra are shown by blue dashed and 
black dotted lines. For convenience the observed spectra, Gaussians, and residual are shifted. Right panel. The same as in left panel, 
but the strong emission line profiles fitted by the single Gaussian. 



Table 1. Emission-Une fluxes 



Line 


F{X) 


1{X) 


Case B" 


Cloudy 


Line 


F(X) 


l(X) 


Case B° 


Cloudy 




a) Near-UV and optical ran 






7281 He: 


1.6+0.4 


0.4+0.1 




0.6 


3726 [O n] 


11.4+0.5 


40.2+1.7 




34.5 


7320 [O II] 


6.2+0.4 


1.5+0.1 




14 


3729 [O II] 


12.9±0.5 


45.6+1.8 




40.4 


7330 [O ii] 


6.1+0.4 


1.4+0.1 




1.1 


3750 H12 


1.0±0.2 


6.5+1.1 


3.07 


3.6 


7751 [Arm] 


7.0+0.4 


1.3+0.1 




1.2 


3770 HI 1 


1.0±0.3 


6.3+1.3 


4.00 


4.5 


8548 Pa 15 


5.4+0.5 


1.1+0.1 


0.53 


0.5 


3797 HIO 


1.5±0.3 


7.8+1.0 


5.34 


5.8 


8601 Pal4 


6.2+0.4 


1.2+0.1 


0.65 


1.0 


3835 H9 


2.4±0.5 


9.6+1.5 


7.37 


7.8 


8753 Pal2 


6.1+0.7 


1.2+0.1 


1.04 


1.2 


3868 [Ne iii] 


7.1±0.7 


21.0+1.5 




60.6 


8865 Pall 


5.7+0.4 


1.1+0.1 


1.35 


1.5 


3889 He i+H8 


4.9+04 


16.6+1.1 


10.60 


17.2 


9017 PalO 


9.2+1.2 


1.5+0.2 


1.80 


2.0 


3967 [Ne iii] 


6.3±04 


16.4+0.9 




18.3 


9069 [S m] 


188.8+6.7 


24.8+0.7 




18.6 


3970 H7 


4.0±0.2 


13.4+0.9 


16.00 


16.5 




b) NIR range 






4026 He i 


1.5+04 


3.7+0.7 




1.8 


9069 [S III] 


118.0+4.2 


15.0+04 




18.6 


4068 [S II] 


1.3+04 


3.0+0.6 




0.9 


9232 Pa9 


20.7+1.8 


2.8+0.2 


2.49 


2.7 


4101 H(5 


9.8+0.6 


23.6+1.2 


26.10 


26.5 


9532 [S m] 


329.1+5.8 


37.6+0.6 




46.1 


4340 Hy 


29.3+1.0 


50.0+14 


47.10 


47.4 


9549 Pa6 


28.8+0.9 


3.7+0.1 


3.57 


3.7 


4363 [O III] 


5.3+0.5 


8.5+0.6 




7.4 


10052 Pa(5 


40.2+2.1 


4.4+0.2 


5.40 


5.6 


4388 He i 


0.9+0.4 


1.3+0.4 




0.5 


10829 He i 


295.1+5.2 


26.8+0.5 




314 


4471 He: 


2.8+0.3 


3.9+0.4 




3.8 


10941 Pay 


76.8+1.5 


7.0+0.2 


8.77 


9.0 


4711 [Any] 


1.6+0.3 


1.8+0.2 




0.7 


12821 Pa/?* 


125.5+2.7 


9.0+0.2 


15.70 


16.0 


4713 He: 


0.9+0.2 


1.1+0.2 




0.5 


15884 Brl4 


4.8+0.3 


0.4+0.0 


0.31 


0.5 


4861 H/3 


100.0+2.0 


100.0+1.3 


100.00 


100.0 


16114 Brl3 


1.8+0.2 


0.2+0.0 


0.39 


0.5 


4959 [0 III] 


283.7+5.5 


259.0+3.4 




220.8 


16412 Brl2 


14.8+0.5 


0.8+0.0 


0.50 


0.6 


5007 [O III] 


817.1+13. 


719.8+8.2 




664.6 


17006 He I 


6.3+0.4 


0.3+0.0 




0.3 


5016 He I 


2.4+0.4 


2.1+0.3 




2.3 


18179 Br9* 


15.7+0.8 


0.8+0.0 


1.21 


1.3 


5876 He i 


22.6+0.6 


11.7+0.3 




9.7 


18756 Pao- 


617.0+11. 


31.0+0.6 


31.90 


32.9 


6300 [O I] 


6.3+0.5 


2.6+0.2 




1.8 


19451 Br5 


34.7+0.9 


1.8+0.1 


1.73 


1.8 


6312 [Siii] 


3.8+04 


1.5+0.1 




1.8 




c) MIR range 






6363 [O I] 


2.0+0.3 


0.8+0.1 




0.6 


10.5 l//m [Siv] 


47.2+15. 


47.2+15. 




84.2 


6548 [N II] 


12.9+0.6 


4.6+0.2 




2.6 


12.81yum [Neil] 


3.9+1.9 


3.9+1.9 




1.1 


6563 Ha 


807.6+14. 


284.8+3.7 


282.00 


285.2 


15.56yum [Ne iii] 


58.9+6.1 


58.9+6.1 




51.2 


6583 [N II] 


37.6+0.9 


13.1+0.3 




7.7 


18.71;um [S III] 


28.9+3.9 


28.9+3.9 




20.7 


6678 He i 


11.0+0.5 


3.6+0.2 




2.6 












6717 [Sii] 


26.5+0.7 


8.5+0.2 




6.1 


c(m 


1.67H 


t0.02 






6731 [Sii] 


21.2+0.6 


7.1+0.2 




5.0 




134+3 






7065 He i 


12.1+0.5 


3.2+0.1 




4.5 


F(H/3)'' 


3.85H 


t0.04 






7136 [Arm] 


23.3+14 


6.0+0.3 




4.8 


EW(abs)'- 


2.0H 


tO.4 







" hydrogen recombination-line relative intensities bv lHumrner &Storevl ( 119871) . 
' affected by the tellmuc absorption. 
' in A. 

'' in units 10"" erg s"' cm"^. 
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Fig. 7. Velocity excess in emission line profiles of a) H/3 /14861; b) [O iii] /14959; c) Ho- A 6563 and d) [O in] A5QQ1. Observational 
data are shown by solid lines and fits by dotted lines (narrow and broad components of lines in blue and summed fits in red). The 
positions of the blue-shifted v'=400 km s ' and red-shifted v=350 km s ' broad components are shown by straight lines. 



and Tf(S in ) were derived from the empirical relations by 
llzotov et"al] ( l2006ah based on a photoionisation model of H 
n regions. We also derive Tf''{0 ii)=(5211+32) K with di- 
rect methods based on extinction-corrected fluxes of emission 
lines from the [O n]/l(3726-H3729)//l(7320-i-7330) emission-line 
ratio, Tf iS in)=(9488+821) K from the [S iii]^9069//l6312 
emission-line ratio and Tf''{S in)=( 13282+449) K from the 
[Siii]^(9069+9532)M63 12 emission-line ratio. We used r,(0 ii) 
for the calculation of O^, and abundances and T^S iii) for 
the calculation of S^^ and Ar-^ abundances. The electron num- 
ber densities Nf'(0 ii)=223+12 and Nf'(S n)=234+65 were 
obtained from the [O n]/}3726M3729 and [S ii]^6717//l6731 
emission-line ratios, respectively. Therefore, the low-density 
limit holds for the H n regions that exhibit the narrow line com- 
ponents considered here. Then, the element abundances do not 
depend sensitively on Ne- The electron temperatures Tf''{0 m), 
Te{0 n), r,(S m), Tf''{S m) and Tf''{0 n), the electron num- 
ber densities Nf''(0 ii) and Nf''(S n), the ionisation correction 
factors (ICFs), and the ionic and total O, N, Ne, S and Ar abun- 
dances derived from the forbidden emission lines are given in 
Table El 

Our oxygen abundance of HG 031203 12+logO/H = 
8.20±0.03 generally agrees with previo us direct determinations 
of 12 -H log O/H in the range ~ 8.0 - 8 2 dProchaska et alil2004 ; 



Hammer et a l.L 1200^ 
2007; LevesQue et al 



iMargutti et al.', '2007*: IWiersema et al 
[2010a; Han et al,, 2010) (Fig.©. In par- 
ticular, the oxygen abundance 12+logO/H - 8.20+0.03 is con 



sistent with 8.19 estimated by ISoUerman et al" (|2005|) from the 
thorough analysis presented bv iProchaska et alT i 2004t) . On the 
other hand, our value is slightly higher than the 1 2 + log O/H 
= 8.12 + 0.04 obtained by Mai-gutti et all (l2007h (VLT opti- 
cal spectra), 8.14 +0.07 bv lHanet all (l2010h and higher than 
the 12 + log O/H = 8.02 ± 0.15 derived bv IProchaska et all 
(120041) (Magellan/IMA CS optical spectra) and 7.96 derived by 
iLevesque etaP (1201 Oah . 



We also recalc ulated 12 + l og O/H with emission - 
line fluxes give r i bv IMargutti et al J (12007). iHanet al.l ( 1201 Ol) . 
IProchaska et al.l (l2004l) and Levesque et alJ ( 2010a ) and obtained 
8.22+0.03, 8.15+0.01, 8.06+0.08 and 8.02, respectively, which 
are higher by ~ 0.01 - 0.10 dex that the values in the original 
papers. 

Abundance ratios log N/O, log Ne/O, log S/O and log Ar/O 
vs. oxygen abundance 12 + log O/H for HG 031203 (this pa- 
per) and for LCGs from SDSS DR7 (llzotov et al.L 12011 ah are 
shown in Fig. |5]by large red filled stars and small blue filled 
circles, respectively. We also plot the available abundance ra- 
tios from the literature and recalculated values from published 
emission-line fluxes by filled and open symbols, respectively. 
Data for HG 031203 are sho wn by stars of d ifferen t colours. 
Specifically, we label data of IMargutti et al J (l2007l) by fi lled 
and open purple stars (HG 031203). iHammer et air(l2006l) by 
filled and open blac k circles (HG 980425 SN and WR re- 
gions, HG 020903), IWiersema etaP (l2007l) b v filled yeflow 
triangles (HG 060218). ILevesque et al l (l2010al) by open light 
blue s quares (HGs 031203, 030329 and 060218), iHan et al.l 
(1201 Oh by open blue rhombs (HGs 31203,030329, 020903 and 
060218). Data for HG 031203 fror n IProchaska et a J (l2004 (us- 
ing solar relativ e abundances by iGrevesse et al.l ( 1996h or by 
iHolwege^ (|200 ih ) are shown by filled green larg e and small stars, 
respec t ively. Abundance ra ti os obtained b y Prochaska et al.l 
( I2()04h . IMargutti et alj ( l2007h . lHammer et alj ( l2006h and our re- 
determined data obtained from their emission line fluxes are 
connected by straig ht lines. The solar abundance ratios by 
lAsplund et al.l (l2009h are indicated by the large purple open cir- 
cles and the associated error bars. We use only those data where 
r^-sensitive [O iii]/l4363A emission line is measured. 

IWiersema et al.l (l2007h found the enhanced N/O ratio for 
several HGs as compared to that in emission-line galaxies. 
However, our value of log N/O for HG 03 1203 ( Fig. [5^) is lower 
by 0.3 - 0.4 dex than log N/O = -0.74+0.20 bv IWiersema et all 



9 - 



N.G.Guseva et al.: A nearby GRB 031203 host galaxy based on VLT/X-shooter observations 

' - Table 2. Physical conditions and element abundances 



8.5 



o 

an 



7.5 



7 - 




LBG 



-15 



-20 



Fig. 8 . Luminosity-metall icity relation for LCGs from SDSS 
DR7 jlzotov et al.L 1201 lab (small blue circles). The location of 
HG 03 1203 is shown by the large filled red star By filled squares 
are denoted the three most metal-poor BCDs (SBS 0335 - 
052E, SBS 0335-052W and I Zw 18) (iGuseva et al.L |2009|) . 
the intermediate-redshift (z < 1 ) extremely low-m etallicity 
emission-line galaxies studied by iKakazu et al.l (l2007h . and lu- 
minous metal-poor star-forming galaxies at z ~ 0.7 studied by 
IHovos et alJ ( 2005h. The po sition of the Lyman-break galaxies at 
z ~ 3 of lPettini et all (1200 ih is marked as LBG with a dotted-line 
re ctangle in the upper right. The emission-line galaxies studied 
bv lGuseva et al.l (l2009l) are shown by dots (red dots are data from 
our observations and black dots are SDSS data). The best linear- 
likelihood fit to the strongly star-forming galaxies derived by 
llzotov et"an (l201 lah is shown by t he black st raight line. The fit to 
emission-line galaxies studied bv lGuseva^ t al. (200^ is shown 



by the red line. Local (z < 0.2) HGs compiled by Margutti et al 



2007h ar e shown by fille d red cir cles. D ata fromlLevesaue et aL 
(I2010ah . iHammer etall (l2006l) and iHanetalJ (l2010l) are 
shown by green, light blue and purple circles, respectively. 

(A colour version of this figure is available in the online jour- 
nal.) 



( I2OO7I) and log N/O = -0.83+0.20 by lProchaska"etaLl ( l200l . 
but it is by . 4 dex higher than log N/O - -1.54±0 . 08 by 
iMargutti et all (l200l . For HG 060218 IWiersema et all (120071) 
obtained log N/O = -1.0±0.4 (our recalculation gives log N/O - 
-1.19+0.88). As Fig.|5^ shows the larger sample of HGs over- 
laps emission-line galaxies from SDSS DR7. For HG 031203 
and HG 060218 and for some emission-line galaxies from SDSS 
DR7 N/O ratio is higher This may be owing to local enrichment 
of the ISM by GRB progenitors or SN remnants. 

From the observed spectrum of HG 031203 we derived log 
Ne/O = -1.15+0.06, using the [Ne iii]/t3868A emission line, 
which is too low compared to the log Ne/O abundance ratio 
in other emission-line galaxies. This is because the flux ra- 
tio [Ne in]/i3868/.i3967 of 1.13 is very low because the [Ne 
ra]/l3868A line is located in a very noisy part of the spectrum 
which renders its flux uncertain. For comparison, the CLOUDY- 
predicted line ratio for [Ne ra]/(^3868)//(^3967) is -3.32 (or 
F{A3868)/ F(A3967) -2.9 3, if C(H/3)=1.67 and the reddening 
function of ICardelli et al.l ((1989) are used). Therefore, we use 
F(3868) = F(3967) x 2.93 for the Ne abundance determination. 
Then log Ne/O = -0.68 +0.05 (see Table Ht, which is consistent 
with the data for LCGs (llzotov et al.L 1201 lah . 



Property 


Value 


Tf (O m)(4959+5007)/4363, K 


12161+307 


TAG II), K 


12034+285 


T (S III) K 


12098+255 


rf XS 111)9069/6312, K 


9488+821 


rf (S iii)(9069+9532)/63I2, K 


13282+449 


Tf(0 ii)(3726+3729)/(7320+7330)), K 


5211+32 




223+12 


A'f XS 11)671 6/6731, cm^^ 


234+65 


0+/H+, (xlO^^) 


1.64+0.14 


02VH+, (xlO^) 


14.20+1.06 






12+log O/H 


8.20+0.03 


N+/H+, (xlO") 


1.55+0.08 


ICFiN)" 


7.82 


N/H, (xlO") 


12.09+0.78 


Ino N/n 

log I^/KJ 




Ne2+/H+, (xlO^) 


3.15+0.28 


/CFCNe)" 


1.04 


Ne/H, (xlO^) 


3.30+0.31 


log Ne/O 


-0.68+0.05 


S+/H+, (xlO«) 


0.24+0.01 


S-+/H+, (xlO^) 


1.64+0.17 


ICF{Sy 


2.09 


S/H, (xlO") 


3.93+0.36 


log S/0 


-1.60+0.05 


Ar2+/H+, (xlO') 


4.49+0.20 


/CF(Ar)" 


1.24 


Ar/H, (xlO') 


5.58+0.23 


log Ar/O 


-2.57+0.04 


" lonisation correction factor. 



IMargutti et al.l (|2007|) do not describe how they derived the 
sulphur abundance. We suspect that their very low sulphur 
abundance shown in Fig. ^ may have been derived because 
the [S iii]/l6312A emission line was not included for the sul- 
phur abundance determination. To demonstrate this suspicion 
we show by smaller open symbols in Fig. |5]: a sulphur abun- 
dance without the [S m]/l6312A emission line. These determi- 
nations result in low log S/O $ -2.0 and are out of the gen- 
eral relation. Note that our recalculated S/ O abundance ra t ios (in 
gala xies with detected [S m U6312A line (IMargutti et al.l (l2007l) 
and'Proc haska et al.l (|2004 for HG 031203 and iHammer et all 
(2006) for HG 980425) are consistent with the S/O determina- 
tions for other emission-line galaxies. 

Overall, element abundance ratios for a sample of HGs and 
our log N/O = -1.12+0.04, log Ne/O = -0.68+0.05, log S/O = 
-1.60+0.05 and log Ar/O = -2.57+0.04 for HG 031203 agree 
well with the abunda nce ratios obtained for a large sample of 
803 spectra of LCGs (llzotov et all 1201 lah . 

3.3. Extinction and hidden star formation 

That the spectrum of HG 031203 has been obtained simultane- 
ously over the entire optical and near-infrared wavelength ranges 
eliminates uncertainties introduced by the use of different aper- 
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tures in the separate optical and NIR observations. In all previ- 
ous studies of low -metallicity emissio n-line galaxies except for 
that of Mrk 59 in llzotoy et al] (l2009l). II Zw 4 0, Mrk 71, Mrk 
99 6, SBS 0335 -052E in llzotov & ThuanI (1201 ih and PHL 293B 
in llzotov et"aD (l201 Ibl) the NIR spectra were obtained in sep- 
arate JHK observations, and there was no wavelength overlap 
between the optical and NIR spectra. The elimination of these 
adjusting uncertainties permits us to compare directly the ex- 
tinctions derived from the optical and NIR spectra. 

In the second column of Table[T]we show the observed fluxes 
F{A) of the emission lines. The extinction coefficient C(Hy6) and 
the equivalent width EW(abs) of hydrogen lines are derived from 
the decrement of hydrogen Balmer lines in the UVB and VIS 
spectra. The extinction coefficient C(Hj6) = 1 .67 is adopted for 
correction of emission line fluxes in all UVB, VIS, and NIR 
ranges. In the third column of the same table the corrected 
fluxes I{X) are shown. In the fourth column of the table the 
th eoretical re combination fluxes for hydrogen lines calculated 
bv lHummer & Sto rey ( 1987) are shown. The relative intensities 
of hydrogen recombination lines dHummer &StorevllT987l) are 
calculated for an electron temperature Te — 12500 K, an electron 
number density A^^. = 100 cm"-*, and the case B theory. 

Note that the comparison of the observed fluxes of Balmer, 
Paschen, and Brackett hydrogen lines after correction for ex- 
tinction with a single value C(Hj6) - 1.67 and the equivalent 
width EW(abs) of h ydrogen lines of 2. OA and the theoretical re- 
combination values ( iHummer &StorevL [19871) shows an agree- 
ment within 2 - 5% for strongest H(5, Hy, Ha, and Paa lines. 
On the other hand, the differences between extinction-corrected 
and predicted intensities for some lines are high (e.g., for Pa/? 
line) because of the strong telluric absorption. The good agree- 
ment implies that a single C(Hy6) can be used over the whole 
~/i/i3200 - 24000A range to correct line fluxes for extinction. 
That the extinction coefficient C(Hy6) does not increase when 
moving from the optical to the NIR wavelength ranges implies 
that ffie NIR emission lines arise in the regions where extinc- 
tion is not as high and that the NIR emission does not have a 
hidden star formation compared to the optical emission lines. 
This appears to be a general result for low-metallicity emission- 
line galaxies (e.g. IVanzi et al 2000, 2002: Izotov et al., ,20091 
1201 IbtEotov & ThuanLl20Tll) . 

To estimate the intrinsic host-galaxy extinction we used 
the far-IR map by Schlegel et al. (1998), from which we in- 
ferred the reddening value for the Milky Way, Ef^wiB - V) ~ 
1.037 [or C(Hj6)mw = 1-52]. We derived the estimate of the 
intrinsic reddening through comparison of C(Hj8)hg=1-67 de- 
rived from Balmer dec rement in the spectrum and C(Hy6)Mw 
= 1.52 bv lSchlegeleFa l. (1998). Thus, the intrinsic C(Hj6)hg 
= 0.15 [or A{V)(HG) = 0.33]. This is typical fo r emission-line 
dwarf galaxies, but lower ffian C(H;6)hg - 0-59 by Margutti et alJ 
(I2007h an d C(H/3)ug = 046 deduced bv lProchaska et alj ( l2004t) 
Recently iHan et all (1201 Oh obtained C(H/?)hg=0.13 + 0.01. 
iLevesque et al.l (1201 Oah derived the total colour excess in the di- 
rection of ffie galaxy C(H/?) = 1.72. If the C(H/3)mw = 1-52 from 
ISchlegel etalT(ll998h is adopted, C(Hj6)hg is 0.2. Both these de- 
terminations are close to our value of intrinsic extinction. 



3.4. CLOUDY stellar photolonlsatlon modelling of the H n 
region and mid-Infrared emission lines 

We next examine the excitation mechanisms of all emission 
lines (not only hydrogen lines) arising in the H n region of HG 
031203. For this purpose, we have constructed stellar photoion- 



isation model for the H ii region using ffie CLOUDY code (ver- 
sion 07.02.01) of lFerlandetd] (Il998h . We list in Table E] the 
input parameters of this model. The first row gives Q{li), the log- 
arithm of the number of ionising photons per second, calculated 
from the extinction-corrected HyS luminosity. The other rows of 
Table[3]list starburst age with the spectral energy distribution of 
ffie ionising steUar rad iation calculated with Starburst-99 models 
dLeithereret al.Lll999h . and Ne, the electron number density, as- 
sumed to be constant with radius. The fourth row gives /, the fill- 
ing factor The remaining input parameters are the ratios of num- 
ber densities of different species to hydrogen, derived from the 
spectroscopic observations in this paper, except for the carbon 
abundance. There are no strong emission lines of carbon in the 
optical and NIR ranges. We therefore calculated its abundance 
using the mean relation of C/O vs. oxygen abundanc e derived for 
low-metallicity H ii regions bv lGarnett et alJ (Il995h . Parameters 
(2(H), Ne and chemical composition are derived from the spectra 
and thus are kept unchanged. The starburst age and filling fac- 
tor are not constrained by the observational data. Therefore, we 
varied these two parameters to achieve the best agreement be- 
tween observations and model predictions. Their adopted values 
are shown in Table [3] 

The CLOUDY-predicted fluxes of the emission lines are 
shown in the fifth column of Table [T] Comparison of the 
extinction-corrected observed and predicted emission-line fluxes 
in both the optical and NIR ranges shows that in general, the 
agreement between the observations and the CLOUDY predic- 
tions is achieved. This implies that a H ii region model includ- 
ing only stellar photoionisation as ionising source is able to ac- 
count for the observed fluxes both in the optical and NIR ranges. 
No additional excitation mechanism such as shocks from stellar 
winds and supernova remnants is needed. 

We showed above that the observed emission-line fluxes in 
the optical and NIR ranges in HG 031203 can be satisfacto- 
rily accounted for with the same extinction coefficient C(Hy6) 
[or the same extinction A(y)]. In other words, there is no more 
hidden star formation in the NIR range as in the optical range 
in regions traced by emission lines, i.e. in ionised gas regions. 
A question then arises: would we see more hidden star forma- 
tion at longer wavelengths, in the mid-infrared range? To inves- 
tigate this question in the MIR r ange, we used the Spitzer obser- 
vations of IWatson et al.l (l2010l) . The MIR emission line fluxes 
relative to the HfS emission Une flux and CLOUDY predicted 
fluxes are shown in Table [T] No extinction correction was ap- 
plied to the observed MIR data because it is presumably small in 
this wavelength range. Comparison of the MIR observed fluxes 
with those predicted by the pure stellar ionising radiation model 
shows agreement to within a factor of ~2 or better. The disagree- 
ment is the worst for the very weak low-ionisation line [Ne ii] 
/1 12. 8 1 yum, whose flux is predicted to be considerably lower ffian 
the observed one. 

The agreement between the observed fluxes of the MIR 
emission-lines and the predicted ones of CLOUDY model, based 
on the extinction-corrected optical emission lines with C(Hj6) = 
1 .67, implies that there is only little more hidden star formation 
seen in the MIR range compared to that in the optical and NIR 
ranges. This implies that in HG 031203 the MIR emission lines 
emerge in relatively transparent regions, which are also seen in 
the optical and NIR ranges. The same conclusion for five other 
emission-line galaxies were made by llzotov et"al] (l2009l) and 
llzotov & ThuanI (1201 Ih . 
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Table 3. Input parameters for stellar photoionisation CLOUDY 
model 



Parameter 


Value 


log Qiuy 


53.73 


Starburst age, Myr 


1.3 




250 


f 


0.01 


log He/H 


-1.09 


logC/H 


-4.30 


logN/H 


-4.92 


logO/H 


-3.80 


log Ne/H 


^.48 


logS/H 


-5.40 


log Ar/H 


-6.37 


"2(H) is the number of ionisin 
''Filling factor. 



Table 4. Dispersions of residual spectra 









h 

s 


Line 




double profile 


single profile 


4861 H/3 


4852-4870 


0.015 


0.034 


4959 [O m] 


4951-4965 


0.039 


0.097 


5007 [O in] 


4999-5013 


0.084 


0.230 


6563 Ua 


6553-6572 


0.076 


0.296 


9532 [S in] 


9520-9547 


0.062 


0.140 


18756 Pao- 


18730-18780 


0.099 


0.151 



"Wavelength range for the 0",^, determination. 
^Dispersions of residual spectra shown in Fig.|6] 



3.5. Kinematic structure 

The brightest emission lines in the spectrum of HG 031203 sig- 
nificantly deviate from a single Gaussian Une profile (Fig. |6] 
right panel). We rule out instrumental effects as the cause be- 
cause no deviations from single Gaussian profiles were detected 
in night-sky line profiles. To study the kinematic structure of 
HG 031203 we reassembled all strong emission lines, obtain- 
ing narrow and broad components separated by -34 km s"'. 
Decomposition of H/3, [O iii] /14959, /15007, Ha, [S iii] /19532 
and Paff /1 18756 emission line profiles into two Gaussian com- 
ponents is shown in Fig. |6] left panel. The observed spectrum 
and the fit are shown by black solid and red dashed lines, re- 
spectively. Two Gaussian components and residual spectra are 
shown by blue dashed and black dotted lines, respectively. For 
convenience the observed spectra, Gaussians, and residuals are 
shifted along the ordinate axis. Dispersions cri es of residual spec- 
tra for two cases: 1) fit of observed profile by two Gaussians and 
2) fit by single Gaussian, are presented in the Table|4] This table 
and Fig.|6]show that the observed bright emission lines are fitted 
much better by the two Gaussians. 

In Table |5] we present observed fluxes, FWHIVI (in km s ') 
and F(A)/F{H/3) ratios of the narrow and broad components for 
14 strong emission lines. The velocity errors are obtained from 
statistical errors of FWHMs from non-flux-caUbrated spectra us- 
ing the IRAF SPLOT routine. 

The FWHlVIs of the narrow-line components and broad-line 
components in the spectrum of HG 031203 are ~90-130 km 
s"' and ~200-330 km s"', respectively. The respective aver- 
age FWHM values for narrow and broad components are ~1 15 
and ~270 km s"'. Note that the Paff profile is likely modi- 
fied by fringes showing a wavy structure (Fig. |6]F). Therefore, 
the decomposition of the Paa line is not as accurate as that 




Fig. 9. Best-fit model SED to the redshift- and extinction- 
corrected observed spectrum. The contributions from the 
stellar and ionised gas components are shown by green 
and blue lines, respectively. The sum of both stellar 
and ionised gas emission is shown by the red line. 
Evidently, the spectrum in the whole wavelength range 
~/l/l3200-24000A is fitted quite well despite the strong ab- 
sorption features in the NIR caused by telluric lines. 

(A colour version of this figure is available in the online jour- 
nal.) 



for other lines. The two velocity components of strong emis- 
sion lines likely arose in two star-forming regions of the host 
galaxy. HG 031203 is a compact galaxy with a major axis of 
which is comparable to the slit width Therefore, it 
is complicated from its kinematical structure to distinguish be- 
tween regular r otating disc motion and stochastic motion o f H 
II regions (e.g. IChristensen et al.L l2008HThone et al.L l2008h. A 
similar kinematic structure was found bylWier sema et alj ( l2007l) 



from two strong emission lines [O iii] /14959, /15007 in HG 
060218. The authors detected two components separated by ~ 
22 km s ' that were repeated in Na i and Ca ii absorption lines . 
The image of AGN candidate Tol 2240-384 dlzotov et al.Ll2010h 
also reveals the two H ii regions, separated by ~80 km s"', if 
the spectral profiles of strongest emission lines are fitted by two 
Gaussians. 

In Fig. |7] the velocity excesses in the emission line profiles 
of H/3 ^4861, [O III] ^4959, Hff A 6563 and [O iii] /15007 are 
shown. Observed profiles are denoted by solid lines and fits by 
dotted lines. Fits of the narrow and broad components are shown 
in blue and summed fits in red. The positions of the blue and 
red velocity excess, corresponding to velocities of v=^00 km 
s"' and v=-i-350 km s"', respectively, are shown by vertical tick 
marks. The bright Paff line is not included in our analysis of the 
velocity excesses because it could be affected by fringes. This 
effect is clearly seen in the Paa profile in Fig. |6] Some traces 
of night-sky absorption lines at the position of Ha (the observed 
wavelength is 7255A) are seen in the blue part of the spectrum 
(Fig.S. 

We cannot obtain physical conditions and heavy element 
abundances for the narrow and broad components separately ow- 
ing to the weakness of the key line [O iii]/l4363. This precludes 
an accurate separation of narrow and broad components of this 
line. There is also insufficient spectral resolution for the decom- 
position of the [O ii]/i3726-i-3729 emission lines into two pairs of 
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Table 5. Observed fluxes of narrow and broad components of strong emission lines 







Narrow 






Broad 




Whole line 


Line 


F(Ar 


FWHM" 


F(A)/F(m 


FW 


FWHM" 




F(A)/F(m 


4340 Hy 


0.802+0.090 


132+10 


0.38+0.04 


0.471+0.074 


330+42 


0.29+0.05 


0.29+0.01 


4861 H/3 


2.100+0.042 


131+3 


1.00+0.03 


1.600+0.027 


290+4 


1.00+0.03 


1.00+0.02 


4959 [O ni] 


6.330+0.121 


131+2 


3.01+0.09 


4.270+0.074 


290+4 


2.67+0.07 


2.81+0.05 


5007 [O III] 


14.700+0.136 


127+1 


7.00+0.16 


12.000+0.117 


279+2 


7.50+0.18 


8.08+0.14 


5876 He i 


0.493+0.041 


126+7 


0.23+0.02 


0.382+0.028 


315+18 


0.24+0.02 


0.22+0.01 


6563 Ua 


13.100+0.096 


106+1 


6.24+0.16 


14.400+0.077 


280+1 


8.57+0.20 


7.99+0.14 


6678 He i 


0.151+0.024 


111+14 


0.07+0.01 


0.256+0.027 


290+30 


0.16+0.02 


0.11+0.00 


6717 S II 


0.525+0.028 


103+4 


0.25+0.01 


0.505+0.028 


261+10 


0.32+0.02 


0.26+0.01 


6731 Sii 


0.358+0.020 


95+5 


0.17+0.01 


0.481+0.023 


257+11 


0.30+0.02 


0.22+0.01 


7065 He i 


0.288+0.021 


117+7 


0.14+0.01 


0.149+0.028 


190+38 


0.09+0.02 


0.12+0.00 


7136 [Arm] 


0.454+0.035 


107+6 


0.22+0.02 


0.452+0.041 


244+16 


0.28+0.03 


0.23+0.01 


9532 [S III] 


6.153+0.094 


117+1 


2.93+0.08 


6.304+0.097 


273+2 


3.94+0.09 


3.25+0.06 


10829 He i 


3.045+0.050 


91 + 1 


1.45+0.05 


7.504+0.093 


221 + 1 


4.69+0.11 


2.92+0.05 


18756 Paa 


12.936+0.067 


124+1 


6.16+0.15 


9.872+0.073 


259+1 


6.17+0.14 


6.10+0.11 


"Observed flux in units 10 er 


I s"' cm"'^. 













*In km s ' . 



Table 6. Parameters of narrow and broad components of strong Balmer lines 





Narrow 


Broad 


Whole 


Line 


F(A)/F{m 


I(A)/I{m 


F(A)/F(m mimp) 


F{A)IF{m l{A)II{Uf}) 


4340 Hy 


0.38+0.04 


0.56+0.05 


0.29+0.05 0.50+0.06 


0.29+0.01 0.49+0.01 


4861 H/3 


1.00+0.04 


1.00+0.05 


1.00+0.03 1.00+0.06 


1.00+0.02 1.00+0.01 


6563 Ha 


6.24+0.16 


2.86+0.07 


8.57+0.20 2.86+0.05 


7.99+0.14 2.85+0.04 


cim 


1.25d 


:0.03 


1.75+0.03 


1.67+0.02 


EW(H/3)" 


70 


+ 1 


56+1 


134+3 


F(H/3)* 


2.10d 


:0.04 


1.68+0.03 


3.85+0.04 


EWCabs)" 


2.0+ 


1.10 


2.0+0.73 


2.0+0.4 



" in A. 

* in units 10"'* erg s"' cm"^. 



narrow and broad components. An inspection of Table |5] shows 
that the observed flux ratios F{A)/F(iy3) are similar for the broad 
components and the entire lines but differ from those of the nar- 
row components. Therefore, we calculated the extinction coeffi- 
cient C{li/3) of Hy, Hfi, and Ho- for the narrow-component region 
and obtained C(Hy6) = 1.24 (Table|6]l. Assuming the same extinc- 
tion for the region with the broad emission, we obtain a broad 
Ha/up flux ratio of ~ 4 or higher than the recombination ratio 
of ~ 2.8, suggesting either a higher density of the broad emission 
region or a higher extinction compared to that in the narrow-line 
region, or Ry different from 3.2. Assuming the higher extinction 
C(Hy6) - 1.75 for the broad-component region, we obtain robust 
estimates of the Balmer-line fluxes (Table|6]l. 



3.6. Luminosity-metallicity relation 



Since the study by iLequeux et al.l d 19791) it was established 
that lower luminosity and lower stellar mass galaxies also have 
lower metallicities. This conclusion has been confirmed using 
large samples of gal a xies such as those from the SDSS (e.g. 
iTremonti et all l2004t iGuseva et al l l2009h . Nevertheless, some 
local low-metallicity star-forming galaxies with strong bursts 
of star formation (e.g. SBS 0335-052E, SBS 0335-052W) and 
some high-redshift galaxies, such as LBGs, deviate from the 
luminosity-metallicity (L-Z) relation for the bulk of emission- 
line galaxies, which are shifted to lower metallicities and/or to 
higher luminosities jKunfli & OstUnll2000l:lGuseva et al.Ll2009t 
llzotov etal.[l2011ah . 



Adopting Z j-22.32 mag from VL T-HFORS2 observation of 
HG 031203 bv llV[arguttietan(l2007h and C(Hj6) = 1.67 we ob- 
tain Bcorr=17.66 mag. At a distance of 430 Mpc taken from 
the NASA/IPAC Extragalactic Database (NED^J we derive an 
absolute magnitude Mb=-20.50. The adopted distance was ob- 
tained from the radial velocity corrected for Virgo infall with 
a Hubble constant of 73 km s ' Mpc For co mparison, the 
absolu te blue m agnitude Mb of HG 31203 from iKewlev et al] 
(l2007h is -19.3. ISavagUo et al.l (l2009h . who estimated Mb in the 
AB system, found -21.11. We compared the absolute magni- 
tudes Mb for HG 31203 wifli the SPS S g absolute magnitudes 
Mg for LCGs from llzqtov et al. ( 2011ah. This is possible thanks 
to the prescriptions of Papaderos et al. (l2008h . who investigated 
the B - g index and concluded that it is very small and vary- 
ing in the range of ~0. 01-0. 03 mag for different star-formation 
histories of a galaxy. 

Fig. [8] shows the relation between the oxygen abundance 
12 + log O/H and the absolute magnitude Mg in the SDSS g 
band f or the extensive sam ple of emission-line galaxies stud- 
ied by iGuseva eTa D (120091) (red dots for our observations and 
black d ots for SDSS sample) and for LCGs from llzotov et"al] 
(1201 lah (smafl blue circles). The location of HG 031203 with 
absolute magnitude and oxygen abundance from this paper is de- 
noted by a large red star. The g apparent magnitudes for LCGs 
and for emission-line galaxies studied by IGuseva et al.l (l2009h 



^ NASA/IPAC Extragalactic Database (NED) is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under con- 
tract with the National Aeronautics and Space Administration. 
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are taken from the SDSS. The correction for extinction was 
made using extinction coefficients C(H/3) of LCGs derived from 
Balmer decrement in the SDSS spectra. Distances of LCGs are 
calculated from redshifts, obtained from strong emission lines. 
With filled black squares we also show the three most metal- 
poor BCDs (SBS 03 35-052E, SBS 0335-052W and I Zw 18) 
(iGuseva et al the intermediate-redshift (z < 1) extremely 

low-m etallicity emission-line galaxies studied by iKakazu et al.l 
(l2007h and lu minous metal - poor s tar-forming galaxies at z ~ 
0.7 studied by iHovos et al.l OOOSh. The area occupied by the 
Lyman-break galaxies at z ~ 3 of Pettini et"an (1200 ih is denoted 
as LBG and is displayed as a dotted-line rectangle. The best 
lin ear-likelihoo d fit to th e strongly star-forming galaxies derived 
by llzotov et"al1 ( 1201 lah is shown in Fig. [8] by a black straight 
line. T he fit to the emission-line galaxies studied bv lGuseva et alj 
(l2009l) is shown by red line. We also show th e location of HGs 
(z < 0.2) compiled by iMargutti et"al] (l2007h by filled red cir- 
cles (980425, 030329, 031203 and 060218). For HG 030329 
and 060218 we used the met allicity obtained from R23 esti- 
mates (Table 9 in the paper of 'Margutti et al., 2007) for which 
the metallicity calibration of IKewlev & Ellison (i20081) was ap- 
phed. Data for HGs 031203, 030329 an d 060218 (green circles) 



are taken from 'Levesque et al.l (l2010ah : for 020903 (fight blue 



circles) from Hammer et alfdZOOel), for HGs 990712, 020903 



03032 9, 031203 and 060218 (purple circles) from iHanetaP 
(120101) . Obviously, HG 03 1 203 in the L - M relation is placed in 
the region of LCGs, which seem to form the bridge between low- 
mass low-metallicity BCD galaxies with extremely high star- 
forming activity and high-redshift LBGs (Izotov et al., 2011a). 

Thus, HG 031203 is very similar in oxygen abundance and 
luminosity to LCGs and together with other HGs lies in the re- 
gion that forms the common lu minosity-metallicity relation, ob- 
tained bv llzotov et al.l (1201 lal) for low-metall i city g alaxies with 
strong star-formation activity. IWatson et all (l201(3l) have con- 
cluded that local BCDs (in meaning, low-metallicity galaxies 
with strong star-formation activity) may be considered as more 
reliable analogues of st ar-formation in t he early universe than 
typical local starbursts. iLevesque et al.l (1201 Oah also came to 
conclusion that LGRB host galaxies appear to fall below L - Z 
relation for local emission-line galaxies, in the region occupied 
by metal-poor galaxies. 

3.7. Star-formation rate 

One of the most important characteristics of galaxy evolutionary 
status is the star-formation rate SFR. We derive the SFR(Hq') us- 
ing the extinction-corrected lumi nosity L{Ha) of t he Ha emis- 
sion line and the relation given by lKennicuttl(ll998h . 



SFR(Ha) 



7.9 X 10"'*'L(Hq'). 



(2) 



In the equation, the SFR is in units of Mq yr L(Hq') 
in erg s corrected for extinction with C(HjS) - 1.67 at a 
distance of 430 Mpc taken from the NED. This results in an 
extinction-corrected Ha luminosity of 7.27 x lO'*' erg s"'. The 
star-formation rate for HG 031203 is then 5.74 Mq yr"', which 
is very clos e to other determination s, for instance, 4.8 Mq yr"', 
obtained bv ILevesque et alJ (l2010ah . The star-formation rate for 
LCGs is in the range 0.7 - 60 Mq yr"' with an average value of 
~4 Mq yr '. Thus, there is no significant difference in the SFR 
for HG 031203 and average value for LCGs. Along the same 
line the range of SFRs in LCGs is compa rable to that in the 
intermediate- r edshif t star-forming g alaxies ( Hoyos et"an , 120051 : 
iKakazu et all l2007h and in LBGs dPettini et all 1200 ih . but is 




Wavelength, S 

Fig. 10. Fraction of gaseous emission to total emission vs. wave- 
length for the modelled spectrum (EW(Hy6) = 134A). The three 
jumps seen at ~/l3660A, ~/18200A and ~/114600A are caused 
by the hydrogen Balmer, Paschen, and Brackett discontinuities 
in the ionised gas emission. 



~10-100 times higher than in typica l nearby BCD s. For example, 
the SFR in 1 Zw 18 is 0.1 Mq yr ' (lThuani l2008l) . 

The iy3 luminosity of HG 031203 derived from the X- 
shooter spectrum and corrected for the extinction with C(Hy6) 
- 1.67 is 2.55 X 10"*' erg s"'. Ba sed on the extensive sample 
of LCGs selected from SDSS DR7. llzotov et al J (|20 1 1 ah showed 
that there is a non-redshift-dependent upper limit of L(Hy6) ~ 
2.5 X 10"*^ erg s"', likely d ue to a self-regulatin g mechanism in 



star formation. All LCGs in llzotov et"aD (1201 lah range in L(H/3) 
from 3 X 10^" to 2.5 x lO'*^ gj.g ^-i -j^^s^ ^j^) jjq 031203 

is in the range of Hj8 luminosities of strongly star-forming galax- 
ies. 

We derived the stellar mass of HG 031203 by fitting its 
spectrum with the stellar po pulations of d i fferent ages. This 
procedure was deve loped by iGuseva et al ] (120061 I2007h and 
llzotov et"al] ( 1201 lah . We include the ionised gas continuum 
emission into the fitting procedure because the neglect of this 
emission in actively star-forming galaxies with EW(Hy6) > lOOA 
leads to overestimates of galaxy ste l lar mas ses by a factor of 
several, as was shown in llzotov et al.l (1201 lah . 

Each fit is performed over the whole observed spectral range 
~AA3200 - 24000A. The shape of the SED depends on several 
parameters. Because each SED is the sum of both stellar and 
ionised gas emission, its shape depends on the relative bright- 
ness of stellar and ionised gas emission. In strongly star-forming 
galaxies, the contribution of the ionised gas can be very high. 
However, the EWs of hydrogen emission lines never attain the 
theoretical values for pure ionised gas emission. This implies 
a non-negligible contribution of stellar emission. We therefore 
parameterise the relative contribution of gaseous emission to 
the stellar one by the equivalent width EW( Hy6). The gas eous 
continuum emission is calculated following lAlleii ( Il984h and 
includes hydrogen and helium free-bound, free-free, and two- 
photon emission. 

The shape of the spectrum also depends on reddening. We 
have no direct observational constraint for the reddening of 
the stellar component, which can be different in principle from 
C(Hj6) obtained from the measured hydrogen line fluxes. 
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Fig. 11. Specific star formation rate SSFR(Ha') vs. the stel- 
lar mass M, is sh own by filled blue circles for LCGs 
dlzotov et all l2011ah . The location of HG 031203 is marked 
by the large red star. The dashed lines correspond from left 
to right to loci with SFR - 1, 10 a nd 10 Mq yr"', re- 
spectively. Data from ISvensson et alJ (1201 Ol) are shown by 
large green open circles. Open gr een and purple stars show 
the position of HG 031203 from ISvensson et al.l ( |2010|) and 
Watson etal. ( 2010|). re spectivelv. Data for different HGs from 
Levesque et al.l ( 2010bh are shown by large filled light-blue 
squares with HG 031203 shown by a light-blue open star). 

(A colour version of this figure is available in the online jour- 
nal.) 



Finally, the SED depends on the star-formation history 
(SFH) of the galaxy. We carried out a series of Monte Carlo sim- 
ulations to reproduce the SED. To calculate the contribution of 
stellar emission to the SE Ds, we adopted a gri d of the Padua 
stellar evolution models by iGirardi et al.l (l2000lfl with a heavy 
element mass fraction Z - 0.0 04. Using these data we calculated 
with the package PEGASE.2 (iFioc & Rocca-Volmeraiig5ll997h 
a grid of instantaneous burst SEDs in a range of ages from 0.5 
Myr to 15 Gyr We adopted a stellar initial mass function with 
a Salpeter slope, an upper mass limit of 100 M© and a lower 
mass limit of 0. 1 M^. Then the SED with any star-formation his- 
tory can be obtained by integrating the instantaneous burst SEDs 
over time with a specified time- varying star-formation rate. We 
approximated the SFH in HG 031203 by a recent short burst, 
which accounts for the young stellar population, and a continu- 
ous star formation responsible for the older stars. The contribu- 
tion of each type of stellar population to the SED is defined by 
the ratio of the masses of the old to young stellar populations, b 
= M(old)/M(young), which we varied between 0.01 and 1000. 

Then the total modelled continuum flux near Hj6 for a mass 
of 1 Mq is scaled to fit the extinction-corrected luminosity of the 
galaxy at the same wavelength. The scaling factor is equal to the 
stellar mass M« of the galaxy. 

The flux ratio of the gaseous continuum to the total contin- 
uum depends on the adopted electron temperature rj.(H^) in the 
H^ zone, because EW(Hj6) for pure gaseous emission decreases 
with increasing r(.(H^). Given the Tg{\\^) is not necessarily 
equal to TJS^ m), we chose to vary it in the range (0.7-1.3)r(,(O 
m). We assume that the extinction coefficient C(Hj6).s,ar.s for the 
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stellar light is the same as C(HyS)g„s for the ionised gas, and vary 
both in the range (0.8-1. 2)C(HySX where C(Hy0) is derived from 
Balmer decrement. 

We ran 5x10^ Monte Carlo models varying simultaneously 
r(young), f(old), h, re(H+), C(H/3)go, and C(H/3)„a„. The best- 
modelled SED is found from minimisation of the deviation 
between the modelled and the observed continuum in ten ranges 
of the spectrum that are free of the emission lines. 

In Fig. |9] we show the best-fit model SED to the redshift- 
and extinction-corrected spectrum. The contributions from the 
stellar and ionised gas components are shown by green and blue 
lines, respectively. The sum of both stellar and ionised gas emis- 
sion is shown by a red line. The spectrum in the whole range 
of wavelengths ~/i/i3200-24000Ais fitted quite well despite the 
presence of the strong telluric absorption features in the NIR part 
of the spectrum. Evidently the contribution of gaseous emission 
is essential because of the EW(H/?) = 134A for the galaxy. 

The fraction of the gaseous emission as a function of wave- 
length in the spectrum of the galaxy is shown in Fig. [10] 
The three jumps at ~/13660A, ~/l8200A and ~/}14600A are 
caused by the Balmer, Paschen, and Brackett discontinuities of 
the ionised gas emission. Between these jumps, the fraction 
of gaseous emission increases with increasing wavelength. As 
Fig.fTOlshows. the fraction of ionised gas continuum-emission in 
the galaxy with EW(H/3) - 134A increases from <10% to -25% 
(in the wavelength range ~4000-8000A), from -14% to -24% 
(in the wavelength range ~9000-15000A) and from -18% to 
-26% (in the wavelength range -15000-20000A). 

The model stellar SED shown in Fig.|9]by the green line is 
used for the stellar mass (M,) determination. The mass of HG 
031203 is M, = 2.5x10^Mq, which is characteristic for a dwarf 
galaxy. Then the specific star-formation rate (SSFR) which is de- 
fined as SSFR(Ha) = SFR(Ha')/M, for the galaxy, is 2.3x10"^ 
yr"'. Ir i Fig. [TT] the SSFR( Ha) vs. the stellar masses M, for 
LCGs (llzotov et al.L 1201 lab are shown by filled blue circles. 
They vary in the range -10^^-10"' yr"'. These values are ex- 
tremely hi gh and similar to t hose found i n high-redshift galaxies 
atz^ 4-6 [stark et al.Ll2009l) and z = 6-8 dSchaerer & de Barrosl 
I2OIOI) . The position of HG 031203 is denoted by a large red 
star The dashed lines correspond from left to right to lo ci with 
SFR - 1, 10 and 100 Mq yr"', respectively. Svensso n et al.l 
(1201 oh obtained the properties (stellar masses, SFRs and SSFRs) 
for 34 HGs, using a large multi-wavelength (0. 45-24^^) dataset 
from GOODS and PANS surveys. Data from ISvensson et al.l 
(I2OIOI) are shown Fig. [TT] by large green open circles. Using 
an unprecedentedly wide wavelength range (-0.3-2xl0''jum), 
IWatson et all ( 120101) have estimated a stellar mass of HG 03 1 203 
as log(M,/Mo) - 9.5. With our estimation of SFR(Ha) we de- 
note the position of HG 031203 in Fig. [TT] (using stellar mass of 
IWatson et al.l (1201 Oh ) by open purple star 

There is an offset between the HG sample and LCGs to 
higher mass (Fig. [TTT i that can be explained by overestimation 
of the galaxy's stellar mass if the ionised gas continuum emis- 
sion is not included in the fitting of the SED. Several HGs have 
lower SFRs (those in the left lower part of Fig.lTTTl. In these cases 
LGRBs arise in ordinary emission-line galaxies with moderate 
star-formation activity. For comparison, the S FR in the pro to- 
type BCD 1 Zw 18 is much lower, 0.1 Mq yr"' (lThuanll2008h . 

It is interesting to note that all parameters of the HG 031203 
host galaxy and other HGs put them into the class of LCGs. It is 
reasonable to adopt that many LGRBs occured in the most lumi- 
nous compact galaxies with low metallicities, because they are 
characterised by the most vigorous star formation. Overall, we 
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find no appreciable difference in heavy element abundances, ele- 
ments ratios, luminosities, star-formation rates and specific star- 
formation rates between HG 031203 and other HGs and LCGs, 
implying their similar nature. 



4. Conclusions 

We have studied the spectrum of the GRB 031203 host galaxy 
(HG 031203) with the VLT/X-shooter spectroscopic observa- 
tions in the wavelength range ~AA32QQ - 24000A . These data 
were co mpared with the data obtained previously bv llzotov et al.l 
(1201 lab for luminous compact emission-line galaxies (LCGs) 
from SDSS DR7. We have arrived at the following conclusions: 
1. We derive the oxygen abundance of 12 + log O/H = 8.20 
+ 0.03 in the H n region of HG 031203. The previou s direct 



determinations dProchaska et al. . 2004t Harnrner et al. . 2006 : 



log 



Margutti et al., '2007; Wie rsema et al r i2oo7n Levesque et al 
2010a; Han et al., 201(3) give the oxygen abundance 12 
O/H in the range ~ 8.0 - 8.2. 

2. We find that the extinction-corrected fluxes of hydrogen 
Balmer, Paschen, and Brackett lines agree well with the theo- 
retical recombination values if a single value of the extinction 
coefficient C(Hy6) = 1.67 is adopted. This implies that there is 
no additional star formation that is seen in the NIR range but is 
hidden in the visible range. The star-forming region observed in 
the optical range is the only source of ionisation. 

3. Using Spitzer MIR emission-line fluxes, we have also 
found that MIR data do not reveal an additional star formation 
that is hidden at shorter wavelengths. Thus, the emission-line 
spectrum of HG 03 1203 in the whole ~ 0.36 - 20 /im wavelength 
range originates in relatively transparent H ii regions. 

4. The profiles of strong emission lines are decomposed into 
two Gaussian narrow and broad components with FWHM ~1 15 
and ~270 km s"', respectively. These components, separated by 
~34 km s likely correspond to two H ii regions with different 
extinction, which is larger in the region with a broad component. 

5. We derive a steUar mass M. = 2.5x10^ Mq for HG 031203 
by fit ting its spect r um with the stellar populations of different 
ages dlzotov et al 1 1201 lat iGuseva et al.l l2006l l2007h . We in- 
clude the ionised gas continuum emission in the fitting procedure 
because the neglect of this emission in the actively star-forming 
galaxy with EW(Hj6 )=134A leads to an overestimate of the M, 
dlzotov et all 1201 lab . 

6. We find that the heavy element abundances, element abun- 
dance ratios, extinction-corrected Ha luminosity L(Ha)=7.27 
X lO"*' erg s"', star-formation rate SFR(Hq') - 5.74 Mq yr"' 
and specific star-formation rate SSFR(Hq') - SFR(Hq')/M, - 
2.3x10"^ yr ' of HG 031203 and other HGs are in the range 
occupie d by the LCGs from SDSS DR7 studied bv llzotov et al.l 
d2011ah In rO m1 /15007/H/? vs. [N ii] i6583/Ha diagnostic di- 
agram dBaldwin et al.L |1981|) and in the luminosity-metallicity 
diagram the GRB host galaxy is also placed in the region of 
the LCGs. This implies that LCGs with extreme star forma- 
tion (containing also green pea galaxies as a subclass of LCGs 
[Cardamo ne et al., 2009i) may predominantly harbour the long- 
duration GRBs. 
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